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ABSTRACT 

We present results for 19 "Lyman Break Analogs" (LBAs) observed with Keck/OSIRIS with an 
AO-assisted spatial resolution of less than 200 pc. We detect satellites/companions, diffuse emission 
and velocity shear, all with high signal-to-noise ratios. These galaxies present remarkably high velocity 
dispersion along the line of sight (~ 70 km s _1 ), much higher than standard star- forming spirals in 
the low-redshift universe. We artificially redshift our data to z ~ 2.2 to allow for a direct comparison 
with observations of high-z LBGs and find striking similarities between both samples. This suggests 
that either similar physical processes are responsible for their observed properties, or, alternatively, 
that it is very difficult to distinguish between different mechanisms operating in the low versus high 
redshift starburst galaxies based on the available data. The comparison between morphologies in the 
UV/optical continuum and our kinemetry analysis often shows that neither is by itself sufficient to 
confirm or completely rule out the contribution from recent merger events. We find a correlation 
between the kinematic properties and stellar mass, in that more massive galaxies show stronger 
evidence for a disk-like structure. This suggests a co-evolutionary process between the stellar mass 
build-up and the formation of morphological and dynamical sub-structure within the galaxy. 
Subject headings: galaxies: kinematics and dynamics - galaxies: starburst - galaxies: evolution 



1. INTRODUCTION 

Our understanding of galaxy formation has changed 
considerably over the course of the last two decades. Re- 
markable progress has been made with numerical simula- 
tions that reproduce the growth of the large-scale struc- 
ture in the universe, and the results from these simula- 
tions agree well with studies of gal axy clusters and the 



Benson & Bower 


2010) 


linear baryonic pJ 


tiysics 



However, the small-scale, non- 



tic structure remains an open question. The so-called 
"gastrophysics" , comprising AGN feedback and super- 
nova winds among other processes, is still poorly under- 
stood. Simulations rely on ad hoc recipes, which are in 
turn based on observational results and are purely phe- 
nomenological; the underlying physical processes are not 
yet known. 

The traditional paradigm of galaxie s forming from 



slowly cooling shock- hea t ed gas (e.g. White & Rees 



Q1978D ; |Mo et al.| fll998| ); |Baughj |2006[ ) and references 



therein) does not seem to apply in many cases. An ele- 
vated fracti on of galaxies at high re dshift display clumpy 
structures ( |Elmegreen et aLl[2008|). which might form 
from internal" instabilities ( |]Noguchi |1999 Immeli et al.| 



2004| Bournaud et al.|2007 ) o r, alternatively, from merg - 
ers of subgalactic gas clumps ( Taniguchi fc Shioya|2001 ), 

~ic 



in agreement with the idea of hierarchical galaxy lorma- 
tion in LCDM models. Furthermore, recent numerical 
simulations indicate that star formation at high redshift 
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might be fed through cooling flows supplying the centers 
of dark matter halo es directly with gas at just below the 
virial temperature flDekel fc Birnboim||2006l |Dekel et al.| 



2009 



Keres et al. 



2009). 



In this context, it becomes 
important to analyze the kinematics in these galaxies, 
and to confront the relative contributions from ordered 
rotation, random motions and merger-induced features 
with predictions from the aforementioned models. Be- 
cause stellar kinematics at high redshift are largely be- 
yond reach of current instruments and telescopes, the 
bright nebular emission line gas is often used as a tracer 
for the underlying kinematics. 

In an early attempt to study kinemat ics of star-forming 
galaxies at z ~ 2 — 3, jErb et al. ( 2006b[ ) analyzed l ong-slit 
spectra of Ha em ission in UV-selected galaxies ( |Steidel| 
et al. 12003, 2004), detecting significant velocity shears m 



12% of the objects in their sample. In all cases, velocity 
dispersion in the ionized gas was high in comparison with 
the observed velocity shears, with v c /a ~ 1. These ob- 
servations are challenging, since they are seeing-limited 
and slit-alignment pl ays an important role in actually 
detecting any sh ears (|Erb et a l. 2006b; Law et al. 2006). 

More recently, |Law et aT ( |2009[ ) improved on this re- 
sult, with spatially resolved kinematics of the gas from 
adaptive-optics (AO) assisted integral- field spectroscopy 
of 12 star- forming galaxies at redshift z ~ 2.5. This tech- 
nique has the advantage of not depending on alignment 
choice, detecting velocity shears all across the extent of 
the galaxy, while the AO system resolves features at sub- 
kpc scales. The authors detect, again, high velocity dis- 
persion values of o~ — 60 — 70 km s _1 . In most cases there 
is no evidence for ordered rotation across the galaxy, and 
in general the gas dynamics appear to be dominated by 
random motions. The authors also find a mild trend 
of rotational properties with stellar mass, with massive 
galaxies typically displaying more pronounced velocity 
shears. 
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In a similar study, Forster Schreiber et al. (2009) stud- 
ied a large sample of 62 star-forming galaxies at simi- 
lar redshifts with the SINFONI instrument. This work 



galaxies, the "supercompact UVLGs", were chosen to 
match those of typical LBGs (e.g., L ~ L^ =3: where 



differs from Law et al. (2009) in that most observations 
are seeing-limited, with spatial resolution elements of ap- 
proximately 4 kpc. In addition, most galaxies in this 



sampl e were drawn from the BzK sample of |Daddi et al. 
( 2004 ) , and are typically two times as massive as the UV- 



selected galaxies. The authors found that their sample 
can be subdivided into three groups: rotation-dominated 
objects, with pronounced velocity shears and v c /o~ val- 
ues of up to 4; dispersion-dominated objects, with little 
to no velocity shear across the major axis; and mergers, 
with multipl e components or peculiar velo city profiles. 
In addition, |Forster Schreiber et al. (2009) also found a 
trend of properties with stellar mass, with more massive 
galaxies presenting higher v c /a ratios and larger sizes. 
A number of observations at intermediate and high red- 
shifts also support the hypothesis of extreme starbursts 
being protodisks resulting either from minor merg ers or 



smooth accretion from the intergalactic medium ( Bouche 



etlI1[20Q7| [Cre sci e t al | [2QQ9| |Wright et al.p)0^ 
pones et al.| (|2010|) also studied the kinematics of the 



ionized gas in high-redshift star-forming galaxies, but 
a sample of strongly lensed objects was used instead. 
The authors were then able to reconstruct the kinematic 
structure by applying models of the gravitational lens, 
achieving much higher spatial resolution (~100 pc) along 
one spatial dimension. Out of a sample of 6 objects, 5 
display characteristics of rotating gas disks, again with 
trends in velocities as a function of size and dynamical 
mass. Although the results help us understand the dy- 
namical structures of such galaxies, it is challenging to 
construct a statistically significant sample of lensed ob- 
jects. Additionally, in many cases the major axis is not 
aligned with the lens shear, in which case the velocity 
shear comprises few resolution elements in the data. 

Studies to date explore complementary regions of pa- 
rameter space. The difference in the prevalence of dif- 
ferent kinematics observed is probably a function pri- 
marily of parent sample, compounded with differences in 
the sensitivity regime of different techniques. The ben- 
efit of AO is that it obtains greater spatial resolution 
but is not sensitive to low surface brightness features (if 
present), while non-AO probes lower surface brightnesses 
and larger radii but with less fidelity. In both cases obser- 
vations are technically challenging, due to the distance to 
the galaxies, which results in low intrinsic spatial resolu- 
tion and cosmic surface brightness dimming. Therefore, 
it is advantageous to observe similar galaxies at lower 
redshifts in order to assess whether certain features de- 
rived from observations at high redshift are intrinsic or 
biased due to observational effects. 

have selected a sample of UV- 

0.2) 



Heckman et al.| 
bright galaxies in t he low-redshift univ erse [z 



from GALEX data (Martin et al. 2005[ ), referred to as 



Ultraviolet Luminous Galaxies (U VLGs). The authors 
found that these galaxies could be subdivided into two 
main groups, one consisting of massive spirals, and the 
other consisti ng of compa ct objects undergoing intense 
starbursts. |Hoopes et al. (2007) further expanded this 
analysis, subdividing UVLGs into three categories with 
respect to their FUV surface brightness. 
The UV characteristics for the most compact of these 



L%_ % is the character istic luminosity of LBGs at z ~ 3). 



Hoopes et al. (2007) found that these objects indeed 
present similar properties to star forming galaxies at 
higher redshift, with comparable star formation rates, 
color s and metallicities, as in ferred from their SDSS spec- 



Basu-Zych et al. (|2007) also determined, from radio 



tra. 

continuum and mid- infrared observations, that these ob- 
jects have significantly less dust attenuation when com- 
pared to galaxies of similar star formation rates in the 
local universe, as is the case for LBGs. We have therefore 
previously referred to these galaxies as "Lyman Break 
Analogs" (LB As), and will do so for the remainder of 
this paper. 

To study the morphologies of LB As, |Overzier et al.| 
20101) obtained HST ultraviolet and optical imag- 



( [20091 



ing of 30 galaxies. In general, their ultraviolet morpholo- 
gies are dominated by clumpy features indicative of mas- 
sive and compact star forming regions, while many fur- 
thermore show clear signs of recent merger events. Inter- 
estingly, when the data is redshifted to z ~ 2 — 4, their 
morphologies are remarkably similar to LBGs at these 
epochs (e.g. [G iavalisco et al. 1996; Papovich et al.||2005 



Lotz et al. 2006; Law et al.|2U07b| ), while the subtle, low 



surface brightness merger features tend to disappear even 
in the deepest rest-frame UV or optical imaging data. 
This implies that on the basis of morphologies alone, it 
cannot be ruled out that LBGs grow through clumpy 
accretion and mergers, perhaps together with rapid gas 
accretion through other means ( jOverzier et~aL]|2QlQ| ). 

Furthermore, strong hydrogen lines and compact sizes 
make them ideal candid ates for IFU spectroscopy. In 
Basu-Zych et al. ( 2009a[ ), we presented preliminary re- 
sults of the IF U survey discussed here for three LB As, 
showing how these galaxies resemble the kinematic struc- 
tures of high-redshift star-forming galaxies. In this work, 
we expand the sample to investigate the ionized gas kine- 
matics of 19 LBAs, observed with spatial resolution down 
to ~200 pc. The paper is divided as follows: in section 
(2), we describe the data acquisition and analysis, includ- 
ing target selection and how we artificially redshift our 
data to z = 2.2 in order to make direct comparisons with 
LBGs; in Section (3), we describe properties of individual 
objects; in Section (4) we describe our results, including 
general trends for these galaxies; in section (5) we discuss 
and analyze the results described in the previous section, 
and in section (6) we summarize our findings. 

Throughout this paper, we assume standard cosmol- 
ogy, with H = 70 km s" 1 Mpc" 1 , ^ m = 0.30 and 
Q A = 0.70. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Sample Selection 

We investigate a subsample of the ultraviolet-luminous 
galaxies (UVLGs). T hese objects were first defined by 
Heckman et al. (2005) to have far-ultraviolet (FUV) lu- 
minosities > 2 x 10 iU L©, which is roughly halfway be- 
tween the characteristic luminosity of present-day galax- 
ies and that of higher redshift Lyman-Break Galaxies 
(LBGs). 

As described in the previous section, |Hoopes et al.| 
( |2007[ ) later expanded the analysis of these objects and 
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subdivided the sample in terms of average FUV surface 
brightness (/1530), using the SDSS w-band half-light ra- 
dius as proxy for the UV size of the galaxies. The sample 
was divided in three categories: large UVLGs (/1530 < 
1O 8 L kpc- 2 ), compact UVLGs (J 153 o > 1O 8 L kpc" 2 ) 
and supercompact UVLGs (/1530 > 10 9 L© kpc -2 ). The 
latter represents the aforementioned LBAs. 

The LBAs are compact systems undergoing intense 
star formation; in fact, they are among the most star- 
forming galaxies in the low-redshift universe. The ob- 
served physical properties, such as metallicity, dust at- 
tenuation, UV/optical morphologies and star formation 
rates, are remarkably similar to those of high-redshift 
LBGs. We further discuss the analogy between low- and 
high-redshift objects in subsequent sections. 

2.2. Observations and Data Reduction 

LBAs are selected to have high surface brightness val- 
ues, which translate into small physical sizes, ranging 
from 0.4 to 1.9 kpc half-light radii in the ultraviolet 
(Overzier et al. 2010). Together with the high star for - 
mation rates up ^ 100 Moyr -1 ( |Hoopes et aL]|2007| ), 
which translates into extremely bright nebular hydrogen 
emission lines, LBAs are highly suitable targets for adap- 
tive optics (AO) assisted integral field spectrogra ph^ 

We have used OSIRIS in the Keck II telescope (Larkin 
et al.|2006[ ). OSIRIS is an integral field unit (IFU) avail- 



able solely for use with AO. It provides a spectral reso- 
lution of R ~ 3800 and a field of view (FOV) of a few 
arcseconds, depending on the configuration utilized. Its 
design is based on a lenslet array, with variable spatial 
pixel scales (spaxels) depending on the need for better 
PSF sampling or a larger FOV. In good weather con- 
ditions, we are able to achieve near diffraction limited 
resolution, or approximately 70 milli-arcseconds (mas) 
FWHM in angular size. 

We have targeted the Pa-a emission line (rest wave- 
length A = 1875.1 nm), which is expected to be ~ 8 
times fai nter than the H-a line, depen ding on gas tem- 
perature (Osterbrock & Ferland 2006). In all cases this 
is redshifted into the redder halt of the K-band, with ob- 
served wavelength varying between 2055 nm < A Q bs ^ 
2350 nm for the objects in our sample. 

The objects observed for this work were selected from 
the original 30 object s observed with HST presented in 
Overzier et al.| ( |2010[ ). Due to a lack of bright nearby 



guide stars, we have used the Laser Guide Star Adap- 
tive Optics (LGS-AO) system for all objects pres ented 



here (van Dam et al.||2006| |Wizinowich et al.|2006[ ). The 
selection of galaxies for each observing run was based 
purely on availability during a given night, proximity of 
prominent sky lines to the wavelength of the Pa-a line 
at each redshift and lesser impact of space command clo- 
sures (when observers are prevented from using the laser 
due to possible collisions with artificial satellites). There- 
fore, no biases were introduced in the data beyond the 
original LBA selection. 

The properties of individual objects are shown in Table 
[l] along with observing information. The stellar masses 
were taken from the SDSS/DR7 MPA-JHU value-added 
catalog Q These masses were calculated by fitting a large 

1 http://www.mpa-garching.mpg.de/SDSS/DR7/ 



grid of spectral synthesis models from Bruzual & Chariot 
(BC03,2003) to the SDSS u' \g' >' \i' >' photometry. The 
lack of near-infrared data and TP-AGB stars from the 
synthesi s library should intro duce an uncertainty of ~0.3 
dex (see |Overzier et aL]|2009[ ), small enough that our re- 
sults, spanning two orders 01 magnitude in stellar mass, 
are unaffected. The BC03 models assume a Chabrier 
(2003) initial mass function, and were chosen to span 
a large range in star formation histories and ages. Prior 
to the fitting, the magnitudes were corrected for emission 
line flux by assuming that the relative contribution of the 
lines to the broad-band photometry is the same inside the 
fibre as outside. Because most of our objects are not or 
barely resolved in SDSS, this is a reasonable assumption. 
The final mass estimate is taken to be the mean value 
of the mass likelihood distribution constructed from all 
models. A more detailed analysis of the masses, ages and 
star formation histories of LBAs based on rest-frame UV 
to far-IR photometry and resolved emission line spec- 
troscopy is currently under investigation. Star formation 
rates presented here are measured from combined Ra and 
MIPS-24/im da ta; they typically pres ent an uncertainty 



up to 0.3 dex flOverzier et al. 2009). For an in-depth 



discussion of proper ties of LBAs and co mpari son with 



high-z galaxies, s ee Ho opes et al. ( |2007| ) and |Overzier 
et^ ( [2^[2010l ). 

Given the limited physical size of the detector, there 
is a trade-off between spatial coverage and wavelength 
coverage; since we are interested in a single emission line, 
we have chosen to use the narrowband mode for most 
galaxies in order to maximize the spatial coverage of the 
data. In most cases we observed with the 50 mas spaxel 
scale; the UV sizes of the remaining objects were larger 
and we chose to use the 100 mas scale with double the 
FOV. 

In many cases, the object occupies a significant por- 
tion of the FOV of the instrument. Because appropriate 
sky subtraction is crucial for a reliable detection of emis- 
sion lines in the data, we have ensured an exclusive sky 
frame was taken in conjunction with each science frame. 
The best strategy to maximize on-target telescope time 
was to observe in 45 minute blocks of science-sky-science 
frames, with 15 minute exposures in each case. Weather 
ranged from acceptable to excellent in all cases, with 
uncorrected seeing (in V band) varying from ~ 1" in 
moderate conditions to 0.5" in the best cases. Weather 
conditions directly affect spatial resolution in our data, 
since the quality of AO corrections depend on the stabil- 
ity and brightness of the laser guide star and the tip-tilt 
star. 

Data were reduced with the OSIRIS pipeline, which 
subtracts the sky frames and translates the two- 
dimensional detector image into a 3D datacube, com- 
posed of two spatial di mensions and one w avelength di- 
mension (for details, see | Wright et al. 12 009 ) . In addition, 
we have written custom 1DL code to further subtract sky 
emission residuals still present in the datacube. This is 
done for each galaxy simply by fitting the ID spectrum 
at spaxels where we believe no signal from the observed 
galaxy exists; this is then subtracted from all spaxels in 
the datacube. 

2.3. Kinematic Maps 
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Table 1 

Summary of LBA Observations 
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a UV half-light radius from HST data 

In order to produce velocity moment maps, we fit 
gaussian functions to the emission lines detected at each 
spaxel. In most cases, our LBA spectra do not show any 
continuum, only the Pa-a line emission. The zero-point 
of the fit is the center of a gaussian fit to the integrated 
one-dimensional spectrum of the collapsed datacube. 

We smooth every datacube spatially with a kernel of 
1.5 — 2 pixels, depending on the data quality and seeing 
in each case. While this results in a slight loss of spatial 
resolution, it also reduces noise, allowing detection of line 
emission at regions with lower surface brightness, espe- 
cially at the outskirts of the galaxies, where gas velocity 
offsets from the center will likely be higher and thus can 
strongly affect our kinematic measurements. In addition, 
to produce the images shown in Fig. [I] we oversample 
the image by a factor of 2, so that features are smoother. 
This is simply a visualization technique and has not been 
used in any of the quantitative analyses discussed in the 
following sections. 

The signal-to-noise (S/N) ratios shown are obtained 
by dividing the area of the gaussian fit to the emis- 
sion line in each spaxel by the sum of the noise fluc- 
tuation over the same wavelength range. The noise is 
determined from a region of the sky with no emission 
line detection. We introduce a minimum threshold of 
S/N = 6 for a fit to be deemed acceptable; anything 
smaller is discarded. This minimizes the presence of ar- 
tifacts in the final maps. This S/N threshold represents 
a detection limit in star formation surface density of or- 
der Ssfr ~ 0.1 M yr -1 kpc~ 2 , comparable to surface 



brigh tness limits determined in |Forster Schreiber et al. 
(2009| and a n order of magnit ude deeper than the data 
presented in |Law et al.| fl2009| ). The velocity-dispersion 
(a) maps, corrected lor instrumental broadening, always 
show values greater than the intrinsic instrumental res- 
olution of ~35 km s _1 , with the exception of some low 
surface brightness spaxels. 

Figure [I] shows the recovered kinematic maps for each 
of the objects in our sample. In each case, the two left 
panels show the HST images of the galaxy, with line emis- 
sion contours overlaid. The third panel shows the zero-th 



moment of the fit, which is simply the total intensity in 
each spaxel, shown as the signal-to-noise of the fit in each 
spaxel. The fourth panel shows the velocity maps, and 
the final panel shows the velocity dispersion maps. We 
also show the resolution element, given by the FWHM of 
a star observed before the galaxy, in the exact same con- 
figuration (band filter and pixel scale). Also shown is a 
horizontal bar indicating a physical size of one kpc at the 
redshift of the galaxy. Three of these galaxies (092600, 
143417 and 210358) hav e been previously analyzed in 
Basu-Zych et all p009a) . 



2.4. Comparison with HST morphologies 

Figure [T] shows the HST images for each galaxy in rest- 
frame optical (left panel) and UV (second-to- left). Im- 
ages are scaled at logarithmic (black) and linear (blue) 
stretch, to distinguish between low surface brightness 
structures and more compact ones. Pa-a flux contours, 
in red, typically enclose 1/3 of the rest-frame optical flux, 
and above 60% of the UV flux. In general we are able 
to detect emission where the bulk of the stellar mass is 
present, unless no significant star formation is present 
(e.g. the southeast components in 080844 and 210358). 

Comparison between both bands in HST shows more 
extended structures in the rest-frame optical, in partic- 
ular at low surface brightness (black). This might in- 
dicate an underlying older stellar population in which 
star- forming regions exist. A complete discussio n of LBA 
morphologies in b oth bands can be found in jOverzier 
eTatn([200m|20T0l). 



2.5. Simulation to High Redshift 

As briefly discussed in section [I] LB As have been de- 
fined on the basis of UV luminosity and surface bright- 
ness thresholds as appropriate for high-redshift Lyman 
break galaxies. Previous studies have supported the 
analogy, finding both apparent and physical properties 
consistent with those of their high-z counterparts. In 
this section we investigate the parallel in terms of gas 
kinematics of LBAs compared to LBGs. 
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Figure 1. We show here the velocity moment maps for all galaxies observed for this work. The two leftmost figures show the HST rest- 
frame optical (left) and UV (right) morphologies, with logarithmic (black) and linear (blue) stretches. The Pa-a S/N levels are overlaid 
in red. There is no UV image available for 101211. The following images show, from left to right, the signal-to-noise ratios, line-of-sight 
velocity in km s _1 and line-of-sight velocity dispersion, also in km s _1 . For the latter two we overplot S/N contours in white. The axes 
show the angular scale in arcsec; orientation is the same in every panel, with north pointing up and east to the left. We indicate in each 
panel the FWHM of a point source as a proxy for spatial resolution and the physical scale corresponding to lkpc at the redshift of each 
galaxy. 
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Figure 1. continued. 



In order to allow for a direct comparison between kine- 
matics of LB A- and LBG-type systems, we have artifi- 
cially redshifted all our galaxies to z ~ 2.^] and reob- 
served the m with the s imula ted IFU prescriptions de- 
scribed by Law et al. (2006). At this redshift, these 



galaxies would be observed m H<x We scale our ob- 
served Pa-a flux maps to the total Ka fluxes determined 
by SDSS. One should notice that the code used to arti- 
ficially redshift our sample represents the exact same in- 
strument, observational setup and reduction software as 



Law et al. (2009[), and has been shown to appropriately 



m 

reproduce actual OSIRIS observations. This ensures the 
robustness of comparisons between th e LBA sample and 
that of LBGs presented in |Law et aT] p009 ) 

We have also artificially redsfiilted our data and sim- 
ulated observations with the SINFONI instrument, in 

2 This precise redshift was chosen to avoid major OH emission 
lines. 




non-AO mode. In this case, the optimal hydrogen line- 
emission surface brig htness detection limits in |Forster 



Schreiber et al.| (|2009) is comparable to our sample: on 
one hand the instrument is more sensitive, Ha is brighter 
and there is no loss due to the adaptive optics system; 
on the other hand, cosmological surface brightness dim- 
ming would make sources up to 200 times fainter per 
solid angle unit. Therefore, we simply degrade our spa- 
tial resolution with a 0.5" gaussian kernel, rebinning our 
datacubes to the nominal 0.125" pixel scale of SINFONI, 
while simultaneously reducing the total angular size of 
the galaxy as determined by the ratio of angular diame- 
ter distances at z = 2.2 and their actual redshift. Exam- 
ples for the resulting velocity maps can be seen in Fig. 

El 



As discussed in Overzier et al. (2010), where a similar 
technique was used lor HST images, much detail is not 
observed due to loss of spatial resolution and/or surface 
brightness dimming. As in the case of HST observations, 
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Figure 1. continued. 

the loss in spatial resolution causes different star-forming 
regions to be confused into one larger clump. This might 
lead to misinterpreting multiple clumps with velocity dif- 
ferences as one larger, smoother rotating disk, with im- 
plications for inferences about its formation mechanism 
(see sections [4JJ5J . This is particularly true for the simu- 
lated SINFONIaata, in which case many LBAs are not 
even spatially resolved. 

These simulations will be used below when comparing 
kinematical measurements of LBAs and actual high red- 
shift galaxies observed. These comparisons, along with 
implications for the analogy between LBAs and star- 
bursts at high redshift, will be discussed in detail in sec- 
tion m 

3. ANALYSIS OF INDIVIDUAL OBJECTS 

In the following sections we briefly describe each ob- 
ject in more detail. Two of these objects (021348 and 
080232) are not resolved even with adaptive optics. They 
present dominant cent ral objects (DCOs) as discussed in 
jOverzier et~aL] Q2009| ). A third object (101211) is too 
faint, and no extended structure is detected. We exclude 
these three objects from the kinematic analysis in subse- 
quent sections. 

3.1. 005527 

This is the only object observed in broadband mode. 
Velocity dispersion is rather uniform across the whole 
galaxy, at about 100 km s _1 . The optical morphology 
is evidently much more extended than the Pa-a emit- 



ting region, which might indicate an underlying, more 
extended, older stellar structure. 

3.2. 015028 

This is an object showing two clearly distinct star form- 
ing regions. There is also a clear velocity shear in the 
east-west direction, which is not aligned with the axis 
connecting the two bright clumps. Velocity dispersion 
is higher in the eastern half of the galaxy. In addition, 
there is some additional emission to the south, at higher 
velocity than the rest of the galaxy; it is unclear whether 
this represents a spiral arm or a tidal tail from an ongoing 
interaction. 

3.3 . 021348 

This is the faintest object observed, and we have only 
been able to detect an unresolved point source in the 
center of the galaxy, in addition to a low S/N region 
(S/N < 10) to the south. It is the first of five objects 
observed with OSI RIS that were class ified as having a 
DCO, according to |Overzier et al.] ( |2009[ ). Since we can- 
not make any inferences about the resolved kinematic 
structure of the galaxy, we have excluded it from any 
further analysis. 

3.4. 032845 

032845 is a bright object, and a significant amount of 
structure is detected. However, velocity shear is remark- 
ably small, and velocity dispersion is, again, relatively 
homogeneous across the galaxy. The HST optical image 
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Figure 2. Velocity maps of 214500 at its intrinsic redshift (left) and artificially redshifted to z ~ 2.2, as observed by OSIRIS (center) and 
SINFONI (right). Legends are the same as in Figure As expected, spatial resolution is lower, and low surface brightness features are 
harder to distinguish. 



shows an antenna-like structure, with distinct nuclei, in 
what appears to be a merger. 

3.5. 035733 

We have been able to detect not only the brightest 
component, but the faint companion to the east, where 
line emission is evidently weaker. A comparison with 
the HST image shows a much more extended structure 
than what is seen here. The western region, however, is 
clearly defined, and shows a definite velocity shear across 
its major axis, resembling a rotating disk, but still with 
line-of-sight velocity dispersion values of approximately 
70 km s _1 , close to the value of the velocity shear across 
the major axis. The companion to the east is at the same 
systemic velocity as the main component. 

3.6. 040208 

This is one of the faintest galaxies we have observed 
(SFR= 2.5M yr _1 ), therefore the signal-to- noise ratio 
is considerably smaller. There are a number of star- 
forming regions northeast of the main component, and 
the velocity offset between them is rather small. 

3.7. 080232 

This is another DCO, like 021348. Again, we detect 
very little emission besides a bright point source in the 
center of the galaxy. This object is also excluded from 
further analysis. 

3.8. 080844 

This is another DCO, but in this case we were able 
to detect emission from the companion to the southeast. 
There is little velocity structure within the main compo- 
nent, but the companion is offset more than 200 km s _1 
from the point source. 

3.9. 082001 

082001 is one of the most elongated objects in our sam- 
ple, which leads to the assumption that it might be disk- 
like structure seen edge-on. The velocity structure seems 
to confirm this hypothesis, with a strong shear across the 
major axis. We are able to detect multiple components, 
indicating there are distinct star-forming regions within 
this disk. 

3.10. 083803 

This object shows a main emission region larger than 
a kpc across, with little velocity structure. In addition, 



we were able to detect emission from a fainter structure 
to the south, with a velocity offset from the main com- 
ponent of ~ 50 km s _1 . This structure is also seen in the 
HST image. 

3.11. 092600 

This is another example of an LBA with a compan- 
ion structure, also evident in the HST image. The com- 
panion presents a ~50 km s _1 shift with respect to the 
main structure. Also evident is a velocity shear across 
the main region itself, albeit small - ~50 km s _1 - es- 
pecially when compared to the velocity dispersion of ap- 
proximately ~100 km s _1 found in the galaxy. This is 
the least massive of our ob jects (log MJMQ = 9.1) and 



has also been described in Basu-Zych et al. (2009a 



3.12. 093813 

This is one of the galaxies with strongest line emission 
in our sample (the Pa-a line is detected at S/N > 400 
in some regions), and therefore we are able to detect the 
substructure with great amount of detail. Multiple com- 
ponents are observed, with velocity offsets greater than 
100 km s _1 between them. Showing signs of a recent 
or ongoing strong merger event in the HST optical data, 
the velocity dispersion seems higher where the merging 
galaxies appear to meet, to the west, where Pa-a emis- 
sion is strongest. 

3.13. 101211 

The emission is weak, and little structure is detected 
beside a faint companion to the northeast. Due to lesser 
data quality in comparison with other galaxies in our 
sample, we do not use this object for our subsequent 
analysis. 

3.14. 113303 

This galaxy shows a remarkable lack of velocity struc- 
ture within the main component, with a shear of a few 
tens of km s _1 , comparable to the instrument resolution 
itself. However, we were able to detect some faint emis- 
sion from a component to the southwest, offset from the 
main region at approximately 100 km s _1 . 

3.15. 135355 

135355 is composed of a large number of small star- 
forming regions, each measuring a few hundred pc across. 
These components show a gradual velocity shear at a 
45 degree angle, indicative of a global velocity structure 
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across the entire galaxy. In addition, there is an elon- 
gated component to the east, visible only in the optical 
HST data and which is likely a merging companion. 

3.16. 143417 

This object presents two clearly distinct regions of star 
formation, along the east-west axis. The regions are at 
distinct velocities with respect to each other. In addition, 
we detect fainter emitting regions to the north and north- 
west, at very different velocities from the two brightest 
regions. These two regions are part of much more elon- 
gated structures, as can be seen in the HST image, which 
shows strong signs of an ongoing interaction. This has 
also been discussed in Basu-Zych et al.| ([2009a ). 



3.17. 210358 

This is the most massive object we have observed, and 
one with very unique f eatures. It is one of the DCO 
objects as described in Overzier et al. (2009), and we 
confirm the existence 01 a bright, unresolved region in 
the center of the galaxy. This region has high Pa-a sur- 
face brightness, with values above 10 -13 erg s _1 cm -2 
arcsec -2 . This galaxy presents the strongest velocity 
shear across its major axis, v s hear ~ 250 km s _1 . This is 



the third object presented in Basu-Zych et al. (2009a). 



3.18. 214500 

This galaxy presents high velocity shear across its ma- 
jor axis, uncommonly so for its low stellar mass (see sec- 
tion 4.2 ). However, its structure is not smooth, and there 
are undetected stellar components to the south, seen in 
the HST image. Likewise, the velocity dispersion map is 
not as well structured as other disk-like galaxies. This 
may indicate a recent merger event. 

3.19. 231812 

This is one of the largest galaxies in our sample, and 
therefore was observed with the 100 mas spaxel scale to 
maximize its field of view. It shows a bright component 
with fainter structure to the south and west. The star- 
forming region to the south has a velocity offset of ~ 75 
km s _1 from the brightest part of the galaxy. 

4. RESULTS 

In this section we discuss some of the analytic results 
obtained from our observations, describing the method- 
ology used to calculate each of the quantities presented. 

4.1. Sizes 

Previous studies of LBAs and high-redshift star- 
bursts infer kpc-scale sizes for the star forming regions, 
from rest-frame UV continuum as observed with HST 



(Overzier et al. 2010, and references therein) and the 
emission l ine regi ons as observed with IFU in struments 
( |Law et al.f2009| |Forster Schreiber et al |2009[ ). Here we 
present our calculations tor sizes of LB As as seen with 
OSIRIS, comparing these figures with results from the 

above-mentioned studies. 

We r eplicate the method described in |Law et al. 



(I2007aj), to allow for a direct comparison with results 
for LBGs utilizing the same instrument. This comprises 
counting the number of spaxels N above a certain S/N 



threshold to represent the size of the star-forming region. 
We use the same threshold, namely S/N > 6. To de- 
termine a radius, we assume galaxies are approximately 

circular, and therefore calculate a radius in spaxels as 

1/2 

r = (N/tt) 1 . This number is corrected for the PSF 



size in each case (see section 2.3) and later converted 
to a physical size at the corresponding spaxel scale and 
redshift of each object. We only use contiguous spaxels 
connected to the brightest region of the galaxy in order 
to exclude companions. Finally, we repeat the process for 
our simulated high-redshift observations of LBAs. Errors 
are typically 0.1 kpc (low-z), 0.4 kpc (OSIRIS) and 2.0 
kpc (SINFONI), given by half the PSF size in each case. 
The results are presented in Table [2] and shown in Fig- 
ure [3] along with actual measurements for high-redshift 
galaxies. Most LBAs present sizes between 1 and 2 kpc, 



consistent with findings from Overzier et al. (2010). 

From Figure [3] we notice this method yields smaller 
sizes at high-redshift in AO-assisted observations; this is 
caused be surface brightness dimming of objects, which 
prevents detection of emission at the outer radii of galax- 
ies. Sizes for our OSIRIS simulated observations are re- 



markably similar to those found in |Law et al. ( 2009 ) for 
LBGs, which have been calculated in an identical man- 
ner; a two-sided Kolmogorov-Smirnov test yields a 97% 
probability of both samples being drawn from the same 
parent population. 

Non-AO observations, however, produce different re- 
sults. Galaxies are apparently larger, in many cases 
due to blending of different components, combined with 
somewhat improved sensitivity at the outer radii, aside 
from obvious loss of resolution. Still, simulated LBAs 



look smaller than ga laxies in the SINS survey ( jForster 
Schreiber et al.||2009[ ) - many, in fact, show sizes smaller 



than the interred uncertainty, which means they are 
essentially unresolved. When comparin g HST sizes of 



LBAs with those found for BzK galaxies, Overzier et al 
( 2010[ ) conclude both samples have similar rest-frame UV 
sizes, but the latter has larger rest-frame optical sizes; 
therefore it is not unreasonable to assume galaxies in 
the SINS survey might be intrinsically larger than both 
LBAs and LBGs. 

4.2. Kinematics and Dynamics of Star Forming 
Galaxies 

The ionized gas in LBAs exhibit very high velocity dis- 
persion^] with median ~ 67 km s _1 and some galaxies 
reaching values above 100 km s _1 . This is much higher 
than those observed in ordinary local star forming galax- 
ies (typical gas velocity dispersions of 5-15 km e.g., 
Dib et al. 2006) but analogous to the increased velocity 
dispersions o bserved in local (ultra-) luminous infrared 
galax ies (e.g. |Arribas et al.||2008| Monreal-Ibero et al.| 
2010). These values are also m good agreement with 
high-redshift star-forming gal axies, as observed bo t h in 



single-slit s pectroscopy (e.g. Pettini et al. 2001| Erb 



l et al. 2006b ) and t he aforementioned integ raTneld studies 
( |Law et al.|2009j |Forster Schreiber et al.|2009[ ). 



3 The global velocity dispersion a measured for the each galaxy 
is an average of each spaxel, weighted by flux. This allows for a 
more accurate measurement than simply measuring the velocity 
dispersion of the whole cube, since it does not incorporate the 
intrinsic velocity shear within the galaxy. 
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Figure 3. Distribution of sizes for LBAs. Black filled histograms 
represent intrinsic values, while red hashed histograms represent 
simulated values at high redshift as would be detected with SIN- 
FONI (top) and OSIRIS (bottom). Blue histograms indicate size 
distributions of actual high-redshift starburst gala xies as measured 
wit h corresponding in struments by Forster Schrei ber et al.| (2009) 
and |Law et aT] ( |2009| >. 

We also measure the velocity shear within each galaxy. 
Since we cannot always precisely define an axis of ro- 
tation, we simply determine the difference between the 
maximum and minimum velocities observed within the 
main body of the galaxy (excluding companions in or- 
der to probe for intrinsic rotation of one star-forming 
region). We determine v max and v m i n as the median of 
the 5-percentile at each end of the velocity distribution, 
so that outliers and artifacts are excluded. The velocity 
shear is then simply defined as i? s hear = ^ ma x — ^min)- 
The values vary between a few tens of km s _1 and over 
200 km s _1 . These measurements are presented in Ta- 
ble [2] In many cases, the velocity shear is not caused 
by actual rotation of the whole galaxy, since there is not 
a significant velocity gradient observed across the entire 
object. 

There is a strong trend of velocity shear with stellar 
mass: more massive objects tend to show greater velocity 
differences between distinct regions of ionized gas. This 
can be seen in detail in Figure Q. Velocity dispersion a, 
also correlates with stellar mass, albeit with a shallower 
slope. For comparison, we also s how in Figure [4] the loca l 
Tully-Fisher relation derived in Bell & de Jong (2001), 
corrected for an avera ge inclination fac tor ol (sin i 



since these objects are not necessarily rotating disks, this 
serves as a comparison with velocity shear in local spirals. 
These values are slightly smaller for a given stellar mass, 
especially at lower masses (up to a factor of 2). Also 
shown is the derived relation fo r star- forming galaxies at 
z ~ 2.2 from |Cresci et al.| ( |2009| ), which shows higher v c - lYC 
values than spirals in the present day; however, in the 
former, the galaxies studied are more massive (M* > 2 — 
3 x 1O 1O M ), and were pre-selected to look like rotating 
disks. 




0.79 (see Appendix in |Law et al. 2009). Although an 



inference for such a relation tor LBAs is not reasonable, 



9.0 9.5 10.0 10.5 n.o 
Log M*/M 

Figure 4. Velocity shear i; s hear {blue circles) and velocity disper- 
sion a (green squares) as a function of stellar mass. The plot shows 
clearly how more massive galaxies show a stronger velocity shear 
than less massive ones, particularly the ones above ~ 10 Mq. 
The same trend, albeit weaker, exists for velocity dispersion a. 
Dashed line shows a power-law fit to our data, wh ile the solid line 
is the Tully-Fisher relation at z <~ according to Bell & de Jong 
(2001). The dotte d-dashed line shows the Tully-Fisher relation at 
z ~ 2 according to |Cresci et al.|(2009| >. 

Due to the difference in slopes, the ratio between veloc- 
ity shear and velocity dispersion (v/cr) is also a function 
of stellar mass (black triangles in Figure [5]). A Spear- 
man's p correlation test shows a ~6% null-hypothesis 
probability of M* and v/a not being correlated. This 
indicates that more massive LBAs have a stronger com- 
ponent of rotational support against gravitational col- 
lapse, as opposed to less massive ones, which are more 
dispersion dominated. 

When artificially redshifted, the v/a ratio decreases, 
from a combination of two effects: on one hand, surface 
brightness dimming causes the high- velocity values at the 
outskirts of the galaxy to be undetected - this is particu- 
larly true for the artificial OSIRIS high-z data (shown as 
red downward triangles in Figure |5l). On the other hand, 
loss of spatial resolution, especially for non-AO obser- 
vations performed with instruments such as SINFONI 
(blue downward triangles in Figure |5| , causes blending 
of features and inner velocity values to dominate, due to 
higher signal-to- noise. The net result is lower v s h ea r val- 
ues. Although our observed v/a values are higher than 
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Table 2 

Kinematic data for LBAs 



Name 


rpa 


V u 

SllGcir 


a 


v/a 


r 

_r ct — OL , IllZ 

(OSIRIS) 


V u 

Sil6£ir , IllZ 

(OSIRIS) 


(J 

(OSIRIS) 


r 

_r ct — OL , IllZ 

(SINF) 


V u 

Sll6£ir , IllZ 

(SINF) 


<7 fai z 

(SINF) 


O (Xyfl 

V W / 


K 

asym 


-^asym,hiz -^asym,hiz 
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1.2 


42 
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0.41 


1.0 


35 


89 


2.5 


18 
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10.2 


0.77 
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1.5 


78 


74 


1.05 
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57 


73 


1.7 


31 


82 


10.0 


0.21 


n i q 
u. iy 
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032845 


1.5 


73 


68 


1.08 


0.6 


13 


46 


3.9 


101 


78 


9.9 


1.60 


0.63 


0.59 


035733 


1.4 


50 


66 


0.76 


1.9 


28 


47 


1.5 


19 


62 


9.9 


0.27 


0.26 


0.25 


040208 


0.6 


53 


50 


1.06 


N/A 


N/A 


N/A 


0.7 


23 


35 


9.3 


0.89 


N/A 


0.44 


080844 


1.2 


27 


92 


0.30 


< 1.5 


16 


95 


3.1 


14 


117 


10.1 


2.16 


0.46 


0.41 


082001 


1.9 


119 


67 


1.78 


1.6 


85 


65 


3.1 


67 


91 


10.0 


0.17 


0.11 


0.08 


083803 


1.4 


41 


49 


0.83 


0.8 


28 


29 


0.9 


13 


45 


9.6 


1.38 


0.50 


0.53 


092600 


1.4 


36 


71 


0.51 


1.3 


23 


54 


1.5 


25 


69 


9.9 


0.61 


0.94 


0.19 


093813 


2.1 


63 


67 


0.94 


1.4 


38 


63 


4.2 


30 


85 


10.0 


0.55 


0.25 


1.85 


113303 


0.7 


14 


30 


0.45 


0.1 


30 


41 


< 2.0 


66 


66 


8.9 


0.66 


0.60 


0.09 


135355 


1.6 


77 


67 


1.15 


1.1 


46 


51 


1.9 


50 


82 


9.9 


0.70 


0.43 


0.19 


143417 


1.7 


73 


67 


1.09 


0.8 


41 


65 


1.9 


37 


67 


10.0 


1.19 


0.21 


0.29 


210358 


1.6 


183 


136 


1.35 


1.0 


72 


161 


3.0 


109 
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10.5 


0.17 


0.12 


0.16 


214500 


1.8 
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55 


1.47 


1.3 


56 
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64 
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9.8 


0.18 


0.20 


0.19 


231812 


2.8 


70 


63 


1.11 


3.6 


55 


55 


2.9 


47 


65 


10.1 


0.28 


0.33 
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high-redshift ones (open circles and squares in Figure J5|, 
when artificially redshifted these galaxies look very sim- 
ilar to high-z star- forming galaxies, with 72% chance of 
being drawn from the same parent population according 
to a standard Kolmogorov-Smirnov test. We present all 
relevant values in Table [2j along with measurements at 
their real redshift. We caution the reader, however, to 
the fact that the observed ratios at low redshift are still 
much smaller than found in local spiral galaxies, which 
have v/a values of 10-20. 

The main kinematic difference when comparing LBAs 
and local spirals comes from gas velocity dispersions, in- 
dicating that LBAs have a dynamically thick structure, 
disk or otherwise. We find it unlikely that the dynam- 
ics in all of the LBAs is actually dominated by rotation, 
given low overall v/a values. Instead, the trend with 
stellar mass might simply indicate a colder, less random 
dynamical structure in the process of forming a disk from 
the dynamically hot gas in more massive galaxies. 

Another quantity one can infer from gas kinematics is 
the dynamical mass of the galaxy, assuming the velocity 
dispersion in the nebular gas is dominated by random 
motions within a gravitational field. In that case, 



Ca 2 r 



(1) 



where G is the gravitational constant, r represents the 
size of the galaxy and C is a proportionality constant 
related to the geometry of the galaxy; for a disk, C = 3.4 , 
whereas for a uniform sphere, C = 5 (|Erb et al.||2006b|). 



Law et al. 



m 



section 



The geometry is not always well determined in LBAs, but 
we assume dispersion-dominated dynamics with C = 5 
to allow for a direct comparison with LBGs as discussed 
(2009). We use Pa-a radius determined in 
ne results are presented in Table [2]). 
Dynamical masses are well correlated witn stellar 
masses. Md yn agrees with stellar masses within a fac- 
tor of two (0.3 dex) in 63% of the galaxies in our sample, 
and they agree within a factor of three for 81% of the 
objects. This means that for most LBAs, the high ob- 
served velocity dispersions can be explained simply by 
random motion of the gas given the observed masses. In 
Figure |6j we present the ratios between dynamical mass 
and stellar mass as a function of stellar mass. We no- 



tice that these ratios are larger for less massive galaxies, 
with the implied dynamical masses an order of magni- 
tude smaller than the observed stellar mass for the most 
massive objects. This supports the hypothesis that the 
dynamical support offered by the rotating disk is more 
significant in the most massive star forming galaxies. It is 
also interesting to notice the same trend, with very sim- 
ilar dynamical and stellar mass values, for high-redshift 
star- forming galaxies, as indicated by red squares. This 
reinforces the analogy between LBAs and LBGs. 

4.3. Kinemetry Measurements 

Another way of assessing the presence of a rotational 
component within the dynamics of the gas in each galaxy 
is provided by the kin emetry method, as introduced by 
Krajnovic et al. (2006). The method comprises a decom- 
position ol the velocity moment maps into its Fourier 
components, that is, for a given ellipse: 

K(i()) = A + Ai sin(^) + B x cos(^) + A 2 sin(2^) + 

+£ 2 cos(2^) + ...,(2) 

where ip is the azimuthal angle along which one measures 
a given velocity moment K (in our case, velocity v or 
velocity dispersion a). Written in another way, 

N 

K(r, = Mr) + Yl k ^ cos t n ^ ~ ^{r))\ , (3) 



where the expansion terms have been redefined as k n 
y/A*TB* and <\> n 



arctan (A n /B n ). For a detailed 



discussion of the met hod, see Krajnovic et al. (2006) and 
Shapiro et all fl2008| . 

For an ideal rotating disk, one would expect the veloc- 
ity profile to be perfectly antisymmetric, that is, the B\ 
term would dominate the Fourier expansion. Likewise, 
the velocity dispersion map is expected to be perfectly 
symmetric, and therefore all terms with the exception of 
Aq would vanish. 



Shapiro et al. (2008) have used this method to ana- 
lyze the dynamics ol high-redshift star-forming galaxies 
observed with the SINFONI instrument. In quantifying 
the asymmetry of the velocity moment maps, they have 
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T LBAs @ z~2 (OSIRIS) 

▼ LBAs @ z~2 (SINFONI) 

O FS+ (2009) 

□ Law+ (2009) 




9.0 9.5 10.0 10.5 11.0 

log M*/M @ 



2.0 



Figure 5. Ratios between velocity shear and velocity dispersion ^shear/^ as a function of stellar mass (black triangles) . The dashed line 
shows a fit to the ^shear/^ data at their intrinsic redshift. We see a moderate trend, indicating more massive galaxies have a stronger 
rotational dynamical component than less massive ones. Also shown as downward triangles are v s h ear /a values for galaxies artificially 
redshifted to z ~ 2 (see text in section [275] ) . In this case, red triangles represent the OSIRIS simulated data, and blue triangles represent 
SINFONI non-AO simulation s. Values irom ac tual high-redshift observations are presented as hollow symbols, representing data from 
Forster Schrei ber et al.| ( |2009| )(open circles) and |Law et al.| ( f2009| > (open squares). The LBAs form an upper envelope with respect to the 
observational data at high redshift, but span a very similar range of parameter space when simulated at z = 2. 
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Log M*/M e 



11.0 



Figure 6. Ratio between dynamical and stellar masses as a func- 
tion of stellar mass. Black filled circles indicate the LBAs in this 
paper, whil e the hollow red squares are the high-z LBGs from Law 
|et al.|j2b09| ). The dashed line indicates M dyri = M*. More massive 
galaxies, both at low and high redshifts, present lower ratios, with 
M dyn ~ 1/10M* for galaxies with logM* ~ lO n M . 

defined the quantities 



^avgjV 



and 



B 



l.v 



(4) 



(5) 



By using local galaxies and numerical models as if ob- 
served at high redshift as templates for disk versus 
merger events, they have found the threshold of 



asym 



asym 



asym 



0.5 



(6) 



to distinguish between rotating disks and mergers. 
Galaxies previously identified by eye as rotating disks 
were correctly classified as disks by the kinemetry 
method, as were galaxies previously identified as merg- 
ers. 



We have used th e same IDL code as presented in Kra- 
jnovic et al. (2006), that, at each semi-major axis, deter- 
mines values tor inclination and ellipticity of the curve 
that will minimize asymmetry. The ellipse center was 
determined as a flux-weighted average of the main body 
of the galaxy, again excluding companions not connected 
to the brightest star-forming region. The ellipses defined 
using the velocity map would then be used with the ve- 
locity dispersion map. 

In Figure [7] we show average kinemetric asymmetry as 
a function of stellar mass. As before, we notice a trend 
in which the most symmetric objects tend to be those 
with high stellar mass. We use the same threshold of 
K asym = 0.5 to distinguish between two categories of 
symmetry. The histogram in the plot shows that more 
asymmetric galaxies (gray bars) are predominantly less 



massive, with one single exception; the symmetric galax- 
ies, on the other hand (green hatched bars), are typi- 
cally more massive, the least massive object having log 
M*/M0 = 9.9. According to a standard Kolmogorov- 
Smirnov test, there is a 0.7% probability that stellar 
masses from K asym > 0.5 are drawn from the same par- 
ent population as the K asym < 0.5 ones. 
Evidently, the threshold of K asym = 0.5 is a simpli 



fication; in fact, an inspection of Figure 7 in |Shapiro 
et al.l p008| shows an overlap of disks and mergers in 



the region where 0.1 < if asym ^ 1-0, which might in- 
dicate instead a transition region between disk galaxies 
and mergers in the kinemetry plot. This region is where 
most star-forming galaxies at high-redshift lie, as is the 
case for the LBAs. 



1.0 



o.i 



i 



3 



-2 







9.0 9.5 10.0 10.5 11.0 
Log Stellar Mass (M@) 

Figure 7. Kinemetric asymmetry measurements as a function 
of galactic stellar mass. Left y-axis shows values of K asym , while 
right y-axis shows quantities for histograms. Gray histogra m shows 
numbe r of g alaxies that would be classified as mergers in JShapiro| 
|et al.| ( |2008t , while the green hatched histogram shows the num- 
ber of galaxies that would be classified as disks. Galaxies with 
high Kasym are predominantly less massive, but the lowest value 
of stellar mass for a galaxy with Kasym < 0.5 is 9.9. 



Finally we compare our resu lts with the high redshift 

In Figure [§1 we show 



simulations from Section 2.5 



the measurements based on the simulations and com- 
pare them to the "intrinsic" values measured in the orig- 
inal (i.e. low redshift) data. The dashed lines show the 
same threshold of K asym = 0.5 used to distinguish be- 
tween disks and mergers. In general, galaxies at high 
redshift present smaller values of if asym, i-e. they appear 
more symmetric than they actually are. One-third of the 
galaxies would be classified differently at high redshift 
(lower-right quadrant). The net effect is that the per- 
centage of galaxies classified as mergers drop from ~70% 
to ~38%. This is a combined effect of signal loss at larger 
radii (where kinematics are less symmetric) and confu- 
sion and blending, smoothing out features that would 
otherwise show departures from a rotating disk. 
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Figure 8. Kinemetry measurements for our high-redshift simula- 
tions as a function of "intrinsic" values measured at low redshift. 



Dashed lines show the same threshold of X a 



: 0.5. Red points 



represent OSIRIS- AO simulations, while blue points represent the 
SINFONI non-AO simulations. The gray shaded area indicates the 
region of the plot where one finds LBAs having high-asymmetry 
values at low redshift but low values a t z ~ 2.2 (lowe r -right quad- 
rant). In the classification scheme of Shapir cTet al.| |2008) these 
objects would likely be classified as rotationally-supported "disks" 



5. DISCUSSION 



In this section we briefly discuss some of the current 
models for galaxy formation, and how they may explain 
the observed properties of both LBAs and LBGs in terms 
of their morphologies and kinematics. In addition, we 
discuss some of the implications of the stellar mass de- 
pendence of observables discussed in the previous section. 

5.1. Ionized gas kinematics as a diagnostic for galaxy 
formation mechanisms 

In light of new techniques and integral-field instru- 
ments, recent studies of the kinematics of ionized gas at 
high redshifts have been used as diagnostics for galaxy 
formation models attempting to explain the distinctive 
properties observed in star-forming galaxies at z ~ 2 — 3. 
In particular, the existence of a large number of rotating 
gas disks with high velocity dispersions at these redshifts 
has been pointed out to support the hypothesis of cold 
gas flows from the IGM directly feeding vigorous star for- 
mation at the center of sufficiently massive dark matter 
haloes (i De kel fc Birnboim||2QQ6| [Dekel et al.|2009| |Keres| 
|et al.|2009 ). The high densities generated would then be 
Toornre unstable, leading to subsequent fragmentation 
into multiple regions and the observed dumpiness of star- 
forming galaxies and, in particular, t he so-called "clump- 



cluster galaxies" at such redsh i fts dlmmeli et al. 



|Elmegreen fc Elmegreen||2QQ5| |Bournaud et al.||2UU8p . 



2004 



These galaxies are ultimately expected to coalesce, with 
individual clumps migrating inwards and creating a bulge 
at th e center of the galaxy (|Noguchil|1999| |lmmeli et al.| 
2004| lElmegreen et al .||2008[ |2UU9[ [Genzel et al.||2UU8[ ). 
Itsnould be noted that the emission lines are produced 



by ionized gas close to star- forming regions, and might 
therefore not be ideal tracers for the dynamics of the 
galaxy as a whole. In fact, comparison between our HST 
rest-frame optical and Pa-a images show we are tracing 
regions that contain approximately a th ird of the total 
stellar mass in the galaxy (section 2.4). Furthermore, 
the gas in these regions is subject to a number of lo- 
cal feedback effects from stellar winds and turbulence, 
and thus may not always repr esent motion of the bulk 
of the dynamical mass within. Lehnert et al. ( 2009| ) ar- 
gue that the high velocity dispersion values could not be 
sustained simply by cosmological gas accretion; instead, 
self-gravity drives the early stages of galaxy evolution un- 
til dense clumps collapse, at which point star formation 
is self-regulated by mechanical output of massive stars. 
Although the inferred dynamical masses from velocity 
dispersion can be explained for the most part as a result 
of rando m m otions within the potential well of the galaxy 
(section 4.2), there is a second-order effect in LBAs that 
supports the influence of star formation-driven turbu- 
lence; in Figure [9] we show v s hear and a as a function of 
star formation rates. Since all variables correlate with 
stellar mass, we present the residuals of a power-law fit 
for all of them with respect to M* , in order to exclude any 
induced correlations. v s hear is independent of star forma- 
tion rates (52% null-hyphotesis probability according to 
Spearman's p test), but more star- forming galaxies show 
stronger velocity dispersion (2% null-hypothesis proba- 
bility), supporting the idea that star formation has an 
effect on a values, generating turbulence from the en- 
ergy output in processes related to star formation. 




-0.8-0.6-0.4-0.2 0.0 0.2 0.4 0.6 

SFR, res 

Figure 9. Velocity shear and velocity dispersion as a function 
of star formation rates. Values are the residuals from a fit with 
respect to stellar mass (see text). 

Another relevant point in this discussion is that rotat- 
ing kinematics do not exclude the possibility of a merger- 
triggered starburst. In a hydrodynamical simu lation of a 
gas-rich major merger, [Robertson et al.|p006[ ) show that 
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for wet mergers rotating disks may form approximately 
~100 Myr after the final coalescence. At this point, a 
large, rotating gaseous disk is formed, with kinemetric 
asymmetry indic es that would in principle rule out the 
merger scenario ( |Robertson fc Bullock|2008| ). The signa- 
ture for interactions might be as subtle as small displace- 
ments of the a map peak from t he center of the galaxy 
( Flores et al.|2QQ6 ). Puech ( 2Q1Q[ ) have used similar argu- 
ments through Lb (J studies ol clumpy, intermediate red- 
shift (z ~ 0.6) galaxies to claim that interactions might 
be responsible for driving star formation at all redshifts. 
In these cases, high gas fractions are a fundamental in- 
gredient in the formation of the disk. Indeed, a num- 
ber of recent CO observations indicate that star-forming 
galaxies at high redshifts have a much large r gas fraction 
than standard spirals in the local universe QTacconi et al. 
20lo] |DaddilFaT1[2QTQl ) 



IFis revealing to compare kinematics of some LBAs in 
our sample with their optical morpho logies as observed 



with HST flOverzier et al.||2009| [20To] ). In FigureJlO 



we 

reproduce the velocity map ol 2T0358 and 1353^5~and 
compare them with their optical image. Although they 
appear as rotating disks in Pa-a, their optical morpholo- 
gies indicate recent major merger events, with quantita- 
tive classification supporting that view. When simulated 
at higher redshift, both the morphology and the veloc- 
ity structure appear much more regular (see Fig. 1 in 



Overzier et al.| ( |2010| ) and Figs. [2] and 10 in this paper). 

Also, on larger scales LBAs seem to support the idea 
of mergers as triggers for the high star formation rates 
observed in these galaxies, as these galaxies tend to pair 
with other galaxies more strongly than a random sample 
does ( |Basu-Zych et al. 2009b). Some studies at high- 
redshilt have concluded that the galaxy number density 
is not high enough to account for mergers in all ob served 

), while 
raction 



starbursts (IConroy et al ||2008| |Genel et al!]|2008 
some morphological studies indicate the merger : 



is high, up to 50%, with M* > 10 *°M^ galaxies under- 
going ~ 4 major merg ers at z > 1 (Conselice et al. 2003; 
|Conselice|[2006l |2009). Furthermore, |Lotz et aEn 2Q08) 
argue that the merger fraction might be even higher, 
since starbursts may outlast morphological asymmetries. 
On the other hand, the aforementioned morphological 
studies are all based on rest-frame UV images; these 
may differ dramatically from the rest-frame optical mor- 
phologies, which better traces the mass distribution of 
the stellar pop ulation - the LBAs themselves presenting 
such contrast ( Overzier et al.||2010[ ). Whether or not the 
same mechanism is triggering star formation at either 
redshift is still unknown. We do not discard the possi- 
bility of an increasing fraction of galaxies at low redshift 
being created by major merger events, and that merg- 
ers of varying mass ratios may be taking place at either 
epoch. 

5.2. The dependence of rotational properties on stellar 

mass 

It has been shown that stellar mass in star-forming 
galaxies at high-redshift correlates with a number of 
physical prope rties, such as metallicity, star formation 
rates and age flErb et al.||2006b|a| jMagdis et al.||2010[ ). 
Some results from kinematic studies ol the ti-a emis- 
sion at high redshift also indicate this dependence, with 
more massive objects being more extended and present- 



ing higher y /a ratios (Law et al. 2009; Forster Schreiber 
|et al.||2009[ ). 

In this work, we have shown that massive galaxies are 
more likely to present disk-like features, as evidenced by 
higher v/a ratios and higher levels of symmetric kine- 
matics, while gas kinematics in less massive objects is 
dominated by random motions, as indicated by higher 
dynamical-to-stellar mass ratios as inferred from veloc- 
ity dispersion measurements. This distinction is particu- 
larly important when taking into account the stellar mass 
function of LBGs at z ~ 2 - 3. |Reddy fc Steidelj ( [20091 ) 
have found that the stellar mass function is particularly 
steep at these redshifts, which means an elevated contri- 
bution from less massive galaxies. That in turn would 
suggest more random dynamics for the majority of star- 
forming galaxies in the early universe, which are respon- 
sible for a significant fraction of the stars observed today 
- ~45% of the present-day stellar mass have formed in 
galaxies with L boi < 10 12 L© ( Reddy fc Steidel||2009| ). 

The dependence on mass is predicted even m more 
traditional star formation models. From a larg e N- 



body/gasdynamical simulation, Sales et al. 



(2009) have 
~ 2 galaxy 



shown that the angular momentum m a z 
depends on halo mass. This dependence extends to stel- 
lar mass in the galaxy, albeit with varying amounts of 
scatter according to feedback efficiency (L. Sales, private 
communication). Whether that can be also a result of 
wet mergers remains unknown. 

It should be noted that the stellar mass presented is the 
global value for the whole galaxy, while Pa-a traces a re- 
gion containing a fraction of the stellar mass. It would be 
interesting to determine whether these relations still hold 
for the stellar mass contained within that small region, 
but for accurate measurements we need high-resolution 
near-infrared imaging, in order to trace stellar mass dis- 
tribution at sub-kpc scales. Alternatively, longer expo- 
sures or more sensitive instruments capable of tracing 
stellar dynamics instead of nebular gas could probe the 
kinematic properties at low surface brightness regions. 
Evidently, this is more difficult at high redshift, where 
cosmological dimming decreases surface brightness val- 
ues by a factor of up to 200. 

It is also unclear whether the trend with stellar mass 
represents an evolutionary effect or simply distinct for- 
mation scenarios. It is tempting to assume these galaxies 
keep forming stars for a period of time, increasing stel- 
lar mass while at the same time settling onto a rotating 
disk. However, this would mean that LBAs would nec- 
essarily keep elevated star formation rates for a period 
over 1 Gyr. Dynamical times of objects containing mul- 
tiple star- forming regions, however, are too short (on the 
order of few tens of Myr), and the galaxy would coalesce 
much more rapidly. Therefore, a continuous inflow of gas 
or a sequence of minor mergers feeding star formation in 
these galaxies would be necessary to keep the observed 
star formation rates. Alternatively, it is possible that 
more massive galaxies have experienced more violent star 
formation episodes in the past, after which the dynami- 
cal structure has cooled down. An in-depth comparison 
with hydro dynamical simulations, with careful examina- 
tion of star formation histories in LBAs, is required to 
examine each hypothesis in detail. 

6. SUMMARY AND CONCLUSIONS 




IFU Spectroscopy of LBAs 



17 




Figure 10. Velocity maps at low (left) and high (center) redshifts for 210358 and 135355. On the right we show the optical morphologies 
of each object as seen by HST, combining optical(orange) and ultraviolet (blue) data. High-z simulated map for 210358 is for OSIRIS data, 
while for 135355 this is the simulated SINFONI data. In the top case we see a galaxy for which a disk is apparent even at low redshift, 
while in the second case we notice the effect of loss of spatial resolution. Both these galaxies are classified as mergers through quantitative 
morphological analysis of the optical images. 



We have performed adaptive-optics assisted observa- 
tions of 19 Lyman Break Analogs (LBAs) with the 
OSIRIS spectrograph at the Keck telescope. By study- 
ing spatially resolved Pa-a emission in these objects, we 
are able to draw the following conclusions: 

(1) All galaxies show high velocity dispersions, indi- 
cating gas dynamics with a strong random component. 
Most galaxies show velocity shears of the same order of 
magnitude as velocity dispersions along the line of sight. 
This is consistent with our general picture of LBAs as 
dynamically young, starburst-dominated galaxies that 
are frequently undergoing mergers as shown by our HST 
data; 

(2) The kinematics in LBAs are remarkably similar to 
high-redshift LBGs, which have also been the target of 
IFU studies. This is demonstrated by artificially red- 
shifting the LBA sample to z ~ 2 and comparing simu- 
lated observations of these galaxies to observations of real 
z rsj 2 L BGs that have been pre sented by other groups 
(see also Basu-Zych et al.||2009a ). All quantitative indi- 
cators of gas properties agree with those found for LBGs 
(e.g. Figures [3] and [5]) . This indicates that our identi- 
fication of LBAs as oeing good local analogs of LBGs 
based on other, previously determined properties (e.g., 
SFR, mass, dust, size, metallicity and morphology) can 
be extended to include their emission line properties as 
well; 

(3) As opposed to IFU observations of high-redshift 
star- forming galaxies, the proximity of LBAs allows for 
a more detailed picture of galactic dynamics. In particu- 
lar, we have high physical resolution and are less subject 
to surface brightness dimming. We show that this bias 
at high can lead to the erroneous classification of sim- 



ilar starburst galaxies at high redshifts, with the kine- 
matic profile appearing smoother and more symmetric 
(e.g. Figure [7]); 

(4) Even in cases where the LBAs resemble a disk at 
low redshift, we cannot rule out mergers as the starburst 
trigger based solely on the gas kinematics. Some disk- 
like galaxies show clear signs of recent interaction in the 
optical imaging data (Figure 10), and the gas disk might 
simply be a result of rapid coalescing of a gas-rich merger. 
Alternatively, a starburst was triggered by a recent infall 
event, while the underlying disk formed previously in a 
more gradual fashion; 

(5) We have shown that whether a galaxy resembles a 
rotating disk depends strongly on stellar mass. The re- 
lationships between stellar mass and disk-like properties 
such as observed velocity shear and dynamical masses as 
inferred from dispersion-dominated gas motions support 
this conclusion. This has strong implications regarding 
the prevalence of disks at high redshift, and might in- 
dicate many of the stars in the local universe have not 
formed in disk-like galaxies. 

Future prospects to study gas assembly in extreme 
starbursts at low and high redshift are excellent. As il- 
lustrated by our work on LBAs, a joint analysis of both 
morphologies and gas kinematics at high resolution and 
sensitivity is absolutely essential for deriving an unam- 
biguous picture of the dynamical state of these systems. 
An important step toward achieving this goal at high 
redshift will be provided by the combination of the ex- 
isting Ha kinematical data with accurate rest-frame op- 
tical morphologies that can now be measured with the 
IR channel on Wide Field Channel 3 aboard HST. Fur- 
thermore, ALMA will allow detection of molecular gas 
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in a large number of star-forming galaxies at redshifts 
z ~ 2 — 3, which should shed more light on the issue 
of gas-rich mergers. At low redshifts, ALMA will allow 
high-resolution measurements of molecular gas distribu- 
tion and kinematics, providing deeper understanding of 
the conversion of gas into stars. Finally, the upcoming 
20- and 30-m class telescopes, which should be opera- 
tional at the end of the decade, will allow IFU studies of 
LBGs with higher sensitivity and resolution levels com- 
parable to what is available now to LB As, while the latter 
will be resolved at scales of giant molecular clouds, and 
we will be able to study the physical processes of star 
formation in situ. In all cases, the LBA sample studied 
in this paper offers a unique low redshift dataset useful 
for contrasting and comparing with starbursts at high 
redshift. 
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